Reversal of spin polarization in Fe/GaAs (001) driven by resonant surface states: 

First-principles calculations 



o 
o 



(N 



CO 



Athanasios N. Chantis, 1 ^ Kirill D. Belashchenko, 2 Darryl L. Smith, 1 

Evgeny Y. Tsymbal, 2 Mark van Schilfgaarde, 3 and Robert C. Albers 1 

1 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
2 Department of Physics and Astronomy and Nebraska Center for Materials and Nanoscience, 
University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA 
3 School of Materials, Arizona State University, Tempe, Arizona 85287, USA 
(Dated: February 1, 2008) 

A minority-spin resonant state at the Fe/GaAs(001) interface is predicted to reverse the spin 
polarization with voltage bias of electrons transmitted across this interface. Using a Green's function 
approach within the local spin density approximation we calculate spin-dependent current in a 
Fe/GaAs/Cu tunnel junction as a function of applied bias voltage. We find a change in sign of 
the spin polarization of tunneling electrons with bias voltage due to the interface minority-spin 
resonance. This result explains recent experimental data on spin injection in Fe/GaAs contacts and 
on tunneling magnetoresistance in Fe/GaAs/Fe magnetic tunnel junctions. 
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Ferromagnetic metal/nonmagnetic semiconductor con- 
tacts have recently attracted significant interest due to 
the possibility to generate non-equilibrium electron spin 
distributions in normal semiconductors and hence be 
practical for spintronics applications The contact 
structures in which electron tunneling dominates trans- 
port properties are used to achieve a sizable spin polar- 
ization of the electric current and produce spin accumu- 
lation in the semiconductor. The spin polarization origi- 
nates from the spin dependence of the wavefunctions and 
densities of states of the ferromagnetic contact. As a re- 
sult, the tunneling transmission coefficients are different 
for majority- and minority-spin electrons. 

Among various ferromagnet / semiconductor structures, 
Fe/GaAs contacts have been extensively studied, showing 
that spin-dependent tunneling through Schottky barriers 
formed by delta doping is an efficient method for inject- 
ing non-equilibrium spin distributions in a semiconduc- 
tor 0, 01 IH . These experiments showed that in biased 
Fe/GaAs contacts the net spin of electrons injected from 
Fe into GaAs is parallel to the magnetization of the fer- 
romagnetic Fe electrode. This implies that majority-spin 
electrons tunnel through the Schottky barrier more effi- 
ciently that minority-spin electrons. 

Recently, however, Crooker et al. || and Lou et al. 
[H observed an anomalous bias dependence of the trans- 
port spin polarization in Fe/GaAs Schottky barrier struc- 
tures. They found that both the magnitude and sign of 
the spin polarization depend on applied bias voltage pro- 
ducing either majority- or minority-spin accumulation in 
GaAs. Moser et al. [7[ observed a related phenomenon 
in Fe/GaAs/Fe magnetic tunnel junctions. They found 
that tunneling magnetoresistance (TMR) changes sign 
with bias voltage, reflecting the reversal of the transport 
spin polarization at the Fe/GaAs interface. To explain 
the experiments by Crooker et al. [5] Dery and Sham Q 



developed a model suggesting that the reversal of spin po- 
larization is due to localized states at the interface formed 
by electrostatic confinement of doping profiles. This ex- 
planation, however, does not take into account a realistic 
electronic structure of the interface which is known to 
be decisive for spin-dependent transport in the tunneling 
regime [9:]. Moreover, it doesn't explain the experiment 
by Moser et al. Q- 

In this Letter, we demonstrate that the observed re- 
versal of the spin polarization in Fe/GaAs(001) tunnel 
contacts is intrinsic to their interface electronic struc- 
ture. The Fe/GaAs(001) interface supports a minority- 
spin interface band lying in the vicinity of the Fermi en- 
ergy Hlj |. This interface band is reminiscent of the 
Fe(OOl) surface band observed experimentally using scan- 
ning tunneling spectroscopy Due to the coupling to 
continuum bulk states in Fe the Fe/GaAs(001) interface 
band evolves into an interface resonant band and strongly 
contributes to the tunneling conductance. The minority- 
spin character of this resonant band leads to the reversal 
of the spin polarization from positive to negative in the 
relevant range of electron energies. This explains the ex- 
perimental findings of anomalous bias dependence of the 
spin polarization in experiments on spin injection Q 
and TMR Q. 

The results of experiments [f| [f| 0] reflect features 
of spin transmission across the Fe/GaAs(001) interface. 
This is due to the transport spin polarization in tun- 
neling geometry being largely controlled by the interface 
atomic and electronic structure . In the case of spin in- 
jection 0,0], electrons injected from Fe into GaAs tunnel 
through the GaAs barrier, then experience scattering by 
a defect or impurity, and further propagate diffusively 
producing spin accumulation in GaAs. Since diffusive 
transport in a nonmagnetic material is independent of 
electron spin, the spin polarization established within 
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GaAs is entirely due to asymmetry in the spin trans- 
mission across the Fe/GaAs(001) interface. A similar ar- 
gument applies to spin extraction from GaAs into Fe. In 
case of magnetic tunnel junctions, variations in TMR are 
expected to reveal spin polarizations of the two interfaces 
However, Moser et al. Q observed a reversed TMR 
only for those Fe/GaAs/Fe tunnel junctions in which one 
interface was " ideal" epitaxial, whereas the other was ei- 
ther oxidized or cleaned by a H + plasma. Therefore, 
their findings reveal features of spin transmission across 
the epitaxial Fe/GaAs(001) interface only. 

To study spin-polarized transport across the 
Fe/GaAs(001) interface we consider a Fe/GaAs/Cu(001) 
tunnel junction with a bec Cu counter-electrode, which 
serves as a detector of spin polarization, in the spirit of 
Ref. [l4j. The bec Cu electrode has a spin-independent 
free-electron-like band structure and a featureless sur- 
face transmission function [13[, making it a perfect spin 
detector. This implies that variations in the spin polar- 
ization of the tunneling current with bias voltage found 
in the calculation performed for the Fe/GaAs/Cu(001) 
tunnel junction are entirely due to the changes in the 
spin transmission across the Fe/GaAs(001) interface. 
This makes the results of our calculations relevant to 
experiments [!, 0, 0] • 

The particular junction studied consists of a semi- 
infinite Fe electrode, 8 monolayers of GaAs barrier, and 
a semi-infinite bec Cu electrode. We consider an As- 
terminated interface, motivated by the experiments on 
spin injection @] where the epitaxial Fe/GaAs inter- 
faces were grown in As-rich environment [la ]. Since in- 
termixing of Fe and As atoms at this interface is not ener- 
getically favorable 11 , 1(| , we assume that the interface 
is abrupt. The small change of the As-Ga interplane dis- 
tance of about 0.14 A due to relaxation [TJJ is not taken 
into account. 

Calculations are performed using the Green's function 
representation of the tight-binding linear muffin-tin or- 
bital (TB-LMTO) method in the atomic sphere approx- 
imation (ASA) 17]. We apply third-order parametriza- 
tion for the Green's function [18J. The electronic struc- 
ture problem is solved within the scalar-relativistic den- 
sity functional theory (DFT) where the exchange and cor- 
relation potential is treated in the local spin-density ap- 
proximation (LSDA). The conductance is calculated us- 
ing the principal-layer Green's function technique 2, 2^ 
within the Landauer-Buttiker approach 2l| . Charge self- 
consistency is achieved before performing the transport 
calculations. 

The spin-dependent transmission coefficient t a (E, k||) 
is calculated for a given spin a =1,1 (where f and I de- 
note majority and minority spin, respectively) as a func- 
tion of energy E and the transverse wave vector k|| which 
is conserved due to the transverse periodicity of the junc- 
tion. The total transmission for a given energy and spin is 
obtained by integrating over kii within a two-dimensional 
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FIG. 1: Results of calculations for a Fe/GaAs/Cu tun- 
nel junction: (a) spin-resolved local DOS for the interface 
Fe monolayer; (b) spin-resolved integrated transmission as a 
function of energy; (c) spin-resolved current density as a func- 
tion of bias voltage; (d) spin polarization as a function of bias 
voltage. The inset shows the spin polarization over an ex- 
tended range of bias [221 ] . In (a) and (b), the Fermi level is 
set at zero energy. 



Brillouin zone (2DBZ): T a (E) = J t a {E^d 2 ^/ (2tt) 2 . 
A uniform 200 x 200 mesh is used for the integration. The 
current density associated with this transmission is ob- 
tained from J a (V) = (e/h) J^ +eV T a (E)dE, where E F 
is the Fermi energy and V is the applied bias voltage. 
This is a reasonable approximation for small voltages 
considered in this work. This definition of J a (V) im- 
plies that for a negative voltage electrons tunnel from 
Fe across GaAs. The spin polarization is defined as 
P=(J T - J X )/(J T + J l ). 

Figs. QJ, and [TJd show the calculated local density 
of states (DOS) at the interface Fe monolayer and the 
integrated transmission as a function of energy for the 
Fe/GaAs/Cu junction. The energies are given with re- 
spect to Ep which is found to be in the middle of the 
GaAs band gap in agreement with previous calculations 



101 111! ]. As is seen from Fig. [TJi, the minority spin 



dominates the interface DOS in the vicinity of the Fermi 
energy throughout the energy interval shown. There is 
a sharp peak in the minority-spin DOS between -50 and 
-160 meV. The majority-spin transmission (Fig. [lb) ex- 
hibits a featureless free-electron-like energy dependence 
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FIG. 3: Ky-resolved minority-spin transmission through a 
Fe/GaAs/Cu (001) tunnel junction for three energies near 
the maximum at E — 106 meV: (a) 45 meV, (b) 106 meV, (c) 
166 meV. 



FIG. 2: Minority-spin Fe local density of states resolved in 
the two-dimensional Brillouin zone by kn with abscissa along 
[10] and ordinate along [01] direction. The upper three panels 
are for E = —121 meV corresponding to the local maximum 
of the minority-spin transmission: (a) bulk, (b) sub-interface 
monolayer, (c) interface monolayer. The lower three panels 
are for E — 106 meV corresponding to the maximum of the 
minority-spin transmission: (d) bulk, (e) sub-interface mono- 
layer, (f) interface monolayer. 

mirroring the featureless majority-spin DOS (Fig. [1^,). 
In contrast, the minority-spin transmission is nonmono- 
tonic and dominates in two energy windows, between 
— 130 and —110 meV and between +50 and +175 meV 
(Fig. [TJd) . The former local maximum corresponds to the 
sharp peak in the minority-spin interface DOS, whereas 
the latter maximum has no distinct analog in the DOS. 

The energy dependence of the transmission is reflected 
in the voltage dependence of the spin-resolved current 
density shown in Fig. [TJ:. It is seen that, while for 
negative bias voltages majority-spin electrons dominate 
the current density, there is a crossover at about +50 
mV which makes the minority-spin current dominating 
at higher voltages up to V = +400 mV (see inset in Fig. 
Q}i). This leads to the reversal of spin polarization at 
about V = +50 mV seen in Fig. p.. At V = +400 mV 
the spin polarization changes sign again reversing from 
anomalous (negative) to normal (positive). The trans- 
mission peak between —130 and —110 meV (Fig. [TJd) 
is too small to change the sign of the spin polarization 
and only leads to a reduction of the spin polarization by 
about 10%. The reversal of the spin polarization with 
bias voltage is the central result of this Letter. In the 
following we will show that an interface resonant band is 
responsible for this anomalous behavior. 

Fig. [2] shows the ky-resolved minority-spin local DOS 
for two monolayers of Fe at the Fe/GaAs (001) interface 
in comparison to the bulk DOS of Fe. The upper three 
panels correspond to the energy E = —121 meV at the 
maximum in the interface DOS (Fig. [1^) and the match- 
ing local peak in the transmission (Fig. Hb). The lower 
three panels correspond to the energy E = +106 meV at 



the maximum in the transmission (Fig. [IJd). It is seen 
that for both energies the interface DOS is strikingly dif- 
ferent from the respective bulk DOS (compare Figs. [5fc 
and 12^,, and Figs. [3J andfSJi). As is evident from Figs. &p 
andf2f, for both energies the interface DOS are character- 
ized by features which have the G% v symmetry intrinsic 
to the atomic structure of the Fe/GaAs(001) interface. 
The topology of these features is preserved at the sub- 
interface monolayer, but their intensity drops down by 
a factor of five (compare Figs. [2t and [2}}, and Figs. [2f 
and [30 • This behavior clearly points to the presence of 
minority-spin interface states at energies E = —121 meV 
and E = +106 meV. The integral DOS for the state at 
E = —121 meV is much higher then that for the state at 
E = +106 meV and consequently the former produces 
the sharp peak in Fig. [T^, whereas the latter is relatively 
broad. 

The analysis of the character of the interface states 
(bands) shows that they arise from a mixture of d 3z 2_ r 2 
and d xy orbitals on the interface Fe sites. These states 
moderately hybridize with bulk Fe minority-spin bands 
and develop into interface resonances. The latter fact is 
evident from their finite width that allows these states to 
be resolved in kn space (see Figs. and[2f). 

The interface resonances contribute to the tunneling 
conductance. The magnitude of their contribution, how- 
ever, strongly depends on their distribution across the 
2DBZ, because the decay of evanescent states in GaAs 
depends on ky . By analyzing the complex band struc- 
ture of GaAs, Mavropoulos et al. [23j demonstrated 
that the decay constant k for the evanescent states has 
a rather deep parabolic global minimum at the T point 
(k|| = 0). This feature strongly suppresses the transmis- 
sion through the resonant states at E = —121 meV, be- 
cause they are located far from the T point (Fig. |5fc). In 
contrast, the resonance at E = +106 meV corresponds to 
the opening of a parabolic pocket at the T point which is 
seen as an ellipse in the surface DOS (Fig. [2f); proximity 
to the f point allows these electrons to tunnel efficiently 
across the GaAs barrier. 

Fig. [3] shows the kii -resolved transmission for three en- 
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ergies around the transmission maximum at E =+106 
meV. It is seen that around the T point the resonance 
band is parabolic and anisotropic, reflecting the Ci v sym- 
metry of the interface. Owing to its location near the 
zone center and a relatively large DOS, this minority-spin 
band dominates the transmission near E=+106 meV. For 
lower energies (Fig. [3^) the resonant band only partially 
crosses the Fermi level (because of its finite linewidth) 
providing fewer states to the tunneling current, while for 
higher energies (Fig. [3);) the crossings occur for larger kn 
which reduces the resonant band contribution due to ky 
filtering in GaAs. This leads to the dominant contribu- 
tion of the minority-spin resonant states in the tunneling 
conductance in a finite energy window (Fig. [T]d) and 
results in the reversal of the spin polarization at bias 
voltages from +50 to +400 mV (Fig. QJi). 

Our results explain the experimental data on spin in- 
jection by Crooker et al. [j| and Lou et al. 6]. In our 
calculations the reversal of the spin polarization occurs 
at positive applied bias voltage. This corresponds to elec- 
trons incoming from GaAs into unoccupied states of Fe, 
that is forward applied bias (spin extraction) in the ex- 
periments [B|, We find positive spin polarization for 
negative applied bias voltage, corresponding to electrons 
incoming from Fe into GaAs, that is reverse applied bias 
(spin injection) in experiments These results are in 

agreement with the experimental data by Lou et al. Q 
(samples A and C). We note, however, that the energy 
of the interface states is sensitive to details of the sample 
preparation and may be affected by inter-diffusion and 
disorder resulting in energy shifts of the order of several 
tenths of an eV. Such a shift may explain why the re- 
versal of the spin polarization occurs when electrons are 
injected from Fe into GaAs in sample B of Lou et al. Q. 

Our results also agree with the TMR data of Moser et 
al. 0- They observed a TMR reversal for Fe/GaAs/Fe 
tunnel junctions with one epitaxial interface and the 
other one either oxidized or cleaned with H + plasma. 
Since no anomalies are observed when both interfaces 
arc disordered, the reversed TMR is entirely due to the 
reversal of the spin polarization at the epitaxial inter- 
face, which occurs at bias voltages from —90 to about 
+400 mV. The minimum in TMR is at V = +50 mV 
which corresponds to electrons transmitted across GaAs 
to the epitaxial interface. This is consistent with our re- 
sults shown in Fig. [TJi. When a sample was annealed 
at 150 °C for lh the first reversal of TMR occured at 
+20 mV instead of —90 mV, while the shape of the trace 
remained essentially unchanged [241 ]. This supports our 
view concerning the sensitivity of interface states to sam- 
ple preparation. 

In conclusion, we have demonstrated that the 
minority-spin resonant states at the Fe/GaAs(001) inter- 
face are responsible for the reversal of spin polarization of 
electrons transmitted across this interface. This explains 
experimental data on spin injection in Fe/GaAs contacts 



and on TMR in Fe/GaAs/Fe magnetic tunnel junctions. 
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